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We theoretically investigated the scheme allowing to avoid destructive space-charge instabilities 
and to obtain a strong gain at microwave and THz frequencies in semiconductor superlattice devices. 
Superlattice is subjected to a microwave field and a generation is achieved at some odd harmonics of 
the pump frequency. Gain arises because of parametric amplification seeded by harmonic generation. 
Negative differential conductance (NDC) is not a necessary condition for the generation. For the 
mode of operation with NDC, a limited space-charge accumulation does not sufficiently reduce the 
gain. 
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There exists a strong demand for miniature, solid- 
state, room temperature operating sources and detectors 
of THz radiation (0.2-10 THz). The need is caused by a 
rapid progress of THz sciences and technologies ranging 
from the astronomy to the biosecurity. Semiconductor 
superlattices (SLs) 0], operating in the miniband trans- 
port regime, are interesting electronic devises demon- 
strating properties of both nonlinear and active media. 
Nonlinearity of voltage-current (UI) characteristic of SL 
gives rise to a generation of harmonics of microwave and 
THz radiation [2J, |3j ; for a high electron mobility and a 
high doping, even such strongly nonlinear phenomena as 
chaos and symmetry breaking || have been predicted. 
On the other hand, Bloch oscillations of electrons within 
a miniband of SL cause an appearance of negative differ- 
ential conductance (NDC) for dc fields (voltages) larger 
than the critical Esaki-Tsu E c Q . Remarkably, the static 
NDC is accompanied by a small-signal gain for ac fields of 
very broad frequency range from zero up to several THz 
0. The feasible device, employing such active media 
properties of SL, is known as the THz Bloch oscillator. 
However, an existence of gain at low frequencies causes 
space-charge instabilities and a formation of high-field 
domains inside SL [f| Q ■ These domains are believed to 
be distractive for the Bloch gain. Recent experiments 
demonstrate that even a use of sophisticated design of 
SL device [si or a development of delicate spectroscopy 
techniques ,9j enable to observe only a very weak Bloch 
gain in dc-biased SLs. 

In these respects important questions arise: Is it pos- 
sible to obtain a gain for a high-frequency field employ- 
ing only a nonlinearity of Ul-characteristic of SL? Can 
we avoid formation of domains in the schemes based on 
nonlinear mechanisms of gain? 

In the present Letter, we examine a feasibility of high- 
frequency superlattice oscillator supporting the regimes 
either with strongly suppressed or even without space- 
charge instabilities. We analyze the situation when a 
pure ac pump field of microwave frequency range is ap- 
plied to SL. Generation arises at third, seventh etc. har- 



monics of the pump. In the weak-probe field limit, the 
mechanism for a gain consists in a combination of para- 
metric amplification of the probe seeded by a frequency 
multiplication of the pump, as well as of nonparametric 
absorption. Remarkably, the oscillator can operate even 
if the amplitude of pump field does not reach E c prevent- 
ing an instability of charge distributions inside SL. Thus, 
we demonstrate the possibility to obtain a gain for a high- 
frequency field employing only nonlinear characteristics 
of SL without using its active medium properties. Effi- 
ciency of the oscillator increases if SL switches partly to 
NDC. For the regime of operation with NDC, we clearly 
show that a space-charge accumulation can only limit but 
not destroy the gain. This mechanism of space-charge in- 
stability stabilization is rather general and applicable to 
other SL devices, like the Bloch oscillator. 

Part of the presented here results has been reported at 
the conferences [To| . 

The paper is organized as follows. (1) Analytic consid- 
eration of a small-signal gain; (2) ID and 3D numerical 
study of a large-signal gain; (3) Effect of a finite Q of res- 
onator on the oscillator efficiency; (4) Analytic analysis of 
space-charge wave stability in the linear approximation 
and within the framework of drift-diffusion model; (5) 
Numerical analysis of nonlinear stage of the instability 
within the drift-diffusion model; (6) Direct comparison of 
the gain found in spatially-homogeneous approximation 
with the gain computed with an account of space-charge 
effects: Calculations using both the ID drift-diffusion 
model and the 3D ensemble Monte Carlo method; (7) 
Effects arising beyond the quasistatic approximation; (8) 
Brief historical review. Some details of the calculations 
are presented in the separate appendixes. 

We consider the response of miniband electrons to the 
action of ac field E(t) — E p + E pr , where E p = E u cos cut 
is the strong pump and the probe, E pr = E n cos uj n t, 
is the nth odd harmonic of the pump, lo„ = nto (n — 
3,5,7,...). SL is placed in a cavity providing a feed- 
back only for the field with frequency u n . In the qua- 
sistatic limit lot < 1, uj n T < 1, the dependence of 
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the drift electron velocity V on the electric field can 
be well described by the Esaki-Tsu formula V(E) = 
2V p (E/E c )/[l + (E/E c ) 2 ], where E c = h/ear is the criti- 
cal field (a is the SL period, r is the intraminiband relax- 
ation time) and V p = V(E C ) is the peak electron velocity 
0. We suppose an operation at room temperature. The 
absorption of the probe field in SL is defined in scaled 
units as A = (v(t) cos(cj„t)), where v = V(E)/V P and 
averaging (...) is performed over the period T — 2tt/uj. 

First, we are interested in a small-signal gain in SL 
(E n < E u ). Expanding v(E) w v(E p ) + v'(E p ) x E pr 
(prime means the derivation with respect to E), we 
present the total absorption A as the sum of Ah = 
(v(E p ) cos(nwt)), A coh — (v'(E p ) cos(2nut)}E n /2 and 
Ainc — {v' (E p ))E n /2 terms, which can be called the har- 
monic, the coherent and the incoherent absorption com- 
ponents, correspondingly (for details see Appendix A). 

We see that while A co h and A inc are dependent on 
both E u and E n , the term Ah is a function of only pump 
strength E Ul and therefore Ah gives the main contribu- 
tion to the absorption of a weak probe. Substituting the 
Esaki-Tsu dependence in the definition of Ah, we find 
A h = {-l) k [2{b-l) 2k+1 ]/bEl k+1 , where 6 = (l+E^) 1 / 2 , 
E^ = E^/Ec and 2k + 1 = n (k = 1, 2,3 . . .). This equa- 
tion describes odd harmonics of the current (c/ Eq. 17 
in 0); in the limit of weak pump E u 1 it takes famil- 
iar form oc E n [2J. Importantly, Ah is negative for the 
odd values of k. Therefore, generated harmonics with 
n = 3, 7, . . . can provide seeding gain for an amplification 
of a probe field. It has no threshold in the amplitude of 
pump E u . Therefore, if the pump amplitude is less than 
the critical field, E^/Ec < 1, the gain at harmonics will 
be not accompanied by the space-charge instabilities in 
SL. 

We turn to the analysis of the term A co h that describes 
a parametric amplification of the probe field due to a co- 
herent interaction of the pump and the probe fields in 
SL. Numerical calculation of the integral demonstrates 
that A co h is always negative for odd n = 3,5,7. In the 
limit of weak pump E u < 1, we find A co h oc E 2n E n . 
Now we can describe the amplification of a weak signal 
at 3a; under the action of a weak pump. Third harmonic 
of a weak pump is generated at cubic nonlinearity of SL 
characteristic, then it can get a seeding gain at the same 
nonlinearity (Ah oc E%), but the next step of paramet- 
ric gain, A co h, uses already the seventh order (~ E 7 ) 
nonlinearity in the y(_E)-dependence. The 5th harmonic 
also can be generated in SL. However, in contrast to the 
case of the 3d harmonic, it cannot be further amplified 
(Ah > 0) and field in the cavity, tuned to 5w, eventually 
evolves to zero. 

Finally, nonparametric effects are described by the 
term A mc = E n /(1 + Elf/ 2 (E n = E n /E c ). We should 
notice that this is the only term which can exist in the 
expression for total absorption of SL in the case u> n /oj is 
not an integer. The absorption Ai nc is always positive. 



We will see that the absorption due to the incoherent 
interactions of fields plays an important role in the sta- 
bilization of space-charge instability in SL. 
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FIG. 1: The dependence of absorption at 3a; on the pump 
amplitude for the different relative probes, n. Inset: Max- 
imal relative amplitude of the probe possessing gain, r/o, as 
function of the pump E^. Calculation with Esaki-Tsu char- 
acteristic (solid) and Monte Carlo technique (symbols). 

To find the efficiency of oscillator we need to consider 
the effects of both a large-signal gain and a finite Q- 
factor of resonator. The stationary value of the cavity 
field, E^, is determined by the balance of gain and loss 
[ill as E*IQf = -A(E U ,E«), where Q' = (u 2 pl T/Lu n )Q, 
E^ = Ef/E c and u pl = (4ire 2 N/em ) 1/2 is the mini- 
band plasma frequency (N is the doping density, mo 
is the electron mass at the bottom of miniband and e 
is the averaged dielectric constant of SL), (for details 
see Appendix B). We found a large-signal absorption 
A(E LU ,E n ) using two methods: (i) simple ID calculations 
with Esaki-Tsu characteristic and (ii) 3D single-particle 
Monte Carlo computations with an account of electron 
scattering at optical and acoustic phonons 01 ■ For the 
Monte Carlo computations we consider GaAs/AlAs SL of 
the period a — 6.22 nm and the miniband width A = 24.4 
meV. Static Ul-characteristic of this SL can be well de- 
scribed by the Esaki-Tsu formula with E c = 4.8 kV/cm, 
V p = 1.44 x 10 6 cm/s and r = 220 fs. We took the 
frequency of pump field lu = 100 GHz. These are our 
default parameters for all computations. 

We start with the case of an ideal cavity (Q — > oo). 
The dependence A(E UI ) for 3o;-generation and for the dif- 
ferent values of relative probe amplitude, r/ = E^/E^, 
are shown in Fig.^ Results of simple ID theory and 3D 
Monte Carlo simulations are in a good agreement. For 
i] = and for small rj <C 1, the dependence of A(E^) fol- 
lows to the corresponding dependencies for the seeding 
gain, Ah(E w ), and the small-signal gain. With a further 
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increase of r), the gain decreases and finally the absorp- 
tion becomes zero for some r/ = i]q: A(E uj , t/q) — 0. The 
value t]q = [E„ /EJ\ 2 determines the maximal SL oscil- 
lator efficiency for the given pump. Inset in Fig. ^ shows 
the dependence of 770 on the pump amplitude. If the am- 
plitude of pump does not reach NDC, E u < E C: the max- 
imal oscillator efficiency is less than 5%. For the oscilla- 
tions in NDC regime, r] 2 reaches 23% for E^jE c w 4. The 
maximal efficiency for 7oj-oscillations is 5% (E^/Ec w 5). 
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FIG. 2: Dependence of the stationary relative amplitude of 
oscillations 770 on Q for the pump amplitudes E^/Ec = 0.95 
(upper) and E^jEc — 5 (lower). Calculation using Esaki-Tsu 
characteristic (solid) and Monte Carlo technique (symbols). 

We also evaluated the oscillator efficiency in the case 
of finite Q. Fig. represents the dependencies of ry = 
Egt/Eu on Q for the doping densities N = 10 15 and 
10 16 cm" 3 . We see that for a heavy doped SL the ef- 
ficiency of generation can reach the maximal efficiency 
(« 23%) even in a cavity with Q < 100 for both opera- 
tional regimes: with NDC and without NDC (c/ Figs. [21 
andQJ. We should also note that oscillator's start up 
also does not require a high-Q resonator because of large 
small-signal gain in SL. 

Although the oscillator can operate even without NDC, 
the maximum of efficiency within spatially-homogeneous 
approximation is reached for E^/Ec > 1. Therefore 
we turn to the consideration of an evolution of space- 
charge instabilities in SL. For lot < 1 a space-time evo- 
lution of the electron density p(x,t), the current den- 
sity j(x,t) and the field E(x,t) inside a SL of the total 
length L, driven by the given voltage U p (t) = U w cos tut 
(U u = E U L), can be well described by the drift-diffusion 
(DD) model 0] . Set of equations consists of the current 
equation j = epV(E) — D(E)dp/dx, the Poisson equa- 
tion dE/dx = 4irer 1 (p — N), the relation to the applied 



voltage U(t) = f E(x, t)dx, and the continuity equa- 
tion edp/dt + dj/dx — [14J. Following the Einstein 
relation the diffusion coefficient is D = kBTV/E. Lin- 
earizing these equations, we find that small fluctuations 
of space-charge with long wavelength will grow, if the 
dielectric relaxation increment, ujd J Q [dV(E p (t))/dE]dt 
(with LJd = w^r), is negative. Taking the integral over 
the period of pump, t = T, we see that the increment 
is proportional to Ai nci which is always positive for all 
n and E^. Therefore the system is stable against small 
fluctuations of charge or field. That is the limited-space- 
charge accumulation (LSA) mode of SL oscillator opera- 
tion. However, in contrast to the traditional LSA mode 
in Gunn diodes driven by dc and ac voltages 15], LSA 
in SL does not require a large amplitude of ac field. 

We also considered the influence of voltage harmon- 
ics calculating numerically the increment for U(t) = 
U p (t) + U n cosuj n t. For the practically important range 
U n < J/*', we found that LSA still works. Next, solving 
numerically DD-model, we consider dynamics of large 
fluctuations. We find that even large fluctuations of 
charge-density are resolved during the time of dielectric 
relaxation oc lu^ 1 . Therefore this type of LSA in SL 
works even in the nonlinear regime. 

Finally, we directly calculated the influence of inhomo- 
geneous distributions of space-charge on a gain in SL. We 
determined the absorption at the nth harmonic of volt- 
age as A d = (((j(x,t)/ j ) cos uj n t)), where j = eV p N 
and averaging ((...)) is performed both over the period 
T and the length L. In the computations a formation of 
domains was caused by a Gaussian spread in the value 
of E c (or r). Domains were periodically created during 
the part of period T when the SL was switched to the 
state with NDC, and then they were annihilated during 
another part of T. We computed the relative decrease 
in gain, 6 = (A — Ad) /A, for different to and U u . De- 
pendence of 5 on Uu/Uc, for a long SL of 130 periods 
and for uj/uid — 0.1, is shown in Fig. [3] Reduction of 
gain due to LSA is 2% for optimal U u jlJ c = 4 and it is 
less than 8% overall. For our default parameters with 
u> = 100 GHz, the value uj/ojd = 0.1 corresponds to the 
doping N = 2 x 10 16 cm" 3 . For the same SL parameters 
but for the pump frequency of tens of GHz, the condi- 
tion ijjjiiJd = 0.1 is satisfied for N ~ 10 15 cm" 3 . As is 
evident from Fig. (lower subplot), resonators with a 
moderate Q still can provide a reasonable efficiency of 
generation in NDC-regime for SLs with the electron den- 
sities N = 10 15 - 10 16 cm" 3 . 

We also found that with an increase of to jtod the value 
of S quickly decreases (see Inset in Fig. |2J). Therefore, 
for higher frequencies or lower doping (ivd oc N) provid- 
ing Lo/ujd > 0.1, the influence of domains on the gain is 
practically negligible. 

We have made additionally the ensemble Monte Carlo 
simulations of space-charge dynamics ^(|. These 3D 
simulations confirmed, in general, our main conclusions 
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FIG. 3: LSA regime of SL oscillator. Relative reduction of 
gain at 3d harmonic 8 vs voltage amplitude U w for io/ujd = 
0.1 and r\ = U^/Uu, = 0.1. Inset: 5 as a function of scaled 
frequency u>/uid for Uu,/U c = 2, r\ = (squares) and r\ — 0.4 
(triangles). 

made in the framework of ID DD-model. Moreover, 
within the Monte Carlo approach we went beyond the 
LSA conditions considering a short (18 periods) but 
heavy doped (N = 10 17 cm -3 ) SL. In this case we found 
that 5 is only w 14% for U u = 4U C . 

In practical terms, the SL oscillator with the paramet- 
ric cascading generation, i.e. the frequency conversion 
followed by the parametric amplification, can operate in 
the whole microwave and THz frequency bands. In the 
present work we restrict our treatment to the microwave 
and low-THz frequency range since there exists plenty of 
serial radiation emitters which can provide a necessary 
power to pump the SL oscillator. Probably the most in- 
teresting and easily realizable case is to use a pump with 
co near 100 GHz in order to obtain a generation in the 
important frequency band of hundreds of GHz ^jj- I n 
order to suppress electrical domains and to use resonator 
with a reasonable Q, the doping of SL should be ~ 10 16 
cm~ 3 . The interaction of miniband electrons with the 
electric field for such frequencies is still quasistatic. 

However, it is instructive to discuss briefly the effects 
arising beyond the quasistatic approximation. Using the 
exact solution of Boltzmann transport e quat ion in the 
constant relaxation time approximation |l8j . we have 
found the following two effects. First, if co is higher than 
some critical co cr , even in an ideal cavity a gain arises only 
if the pump E u is larger than some threshold value Eth- 
For the gain at third harmonic io cr w 0.34r _1 (« 250 
GHz); importantly, we have E t h < E c for the frequen- 
cies satisfying to < 0.38t _1 (« 270 GHz). Therefore, the 
generation of radiation at 3lo without NDC and, as con- 
sequence without an influence of domains, is possible for 
the pump field with frequencies co < 270 GHz. The sec- 



ond effect, which can arise beyond the quasistatic ap- 
proximation, is the possibility to have a gain at the 5th 
and the 9th harmonics. For example, a small gain at 5th 
harmonic arises for to > 0.14r _1 (« 100 GHz). 

Theoretical research devoted to a gain arising at fre- 
quency multiplication in SLs has some history. To the 
best of our knowledge, Pavlovich presented the first cal- 
culations showing a principle possibility of a small-signal 
gain in the presence of a strong pump with lot > 1, i f the 
frequency of the probe is some half-integer of to |l9t |2fJ • 
Romanov pointed out that such a gain should originate 
from a multiphoton parametric amplification (21j. Pos- 
sibility to get a gain at harmonics of microwave pump 
also has been discussed in j^, Because references 
to each other are absent in the papers [lil EH E2T | , we 
suppose that all these works were done independently. 
This interesting activity did not receive much attention 
so far. We should also underline that a possible role of 
space-charge instabilities in SL has not been analyzed in 
these previous contributions. 
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APPENDIX A. WEAK-SIGNAL ABSORPTION 
IN THE QUASISTATIC LIMIT 

Here we derive the expression for absorption of the 
weak probe field E pr = E n cos{ntot) in the presence of 
strong pump field E p = E u cos cot, (E n <C E ). 

Introduce the Fourier expansions for the electron ve- 
locity V(E(t)) and its derivative V'(E(t)) (prime means 
derivation with respect to E) as 

1 °° 

Y(E pump (t)) = -f + ^ cos(lcot), (1) 
i=i 

where Y stands either for V or for V'. Fourier transform 
Yi is defined as 

% = % [ Y(E p (t))cos(lLot)dt (2) 
1 Jo 

Substituting the Fourier expansion JQ) and the Taylor 
expansion 

V(E) « V(E p ) + V'(E p ) x E pr 
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in the definition of absorption, we obtain 



A = Vp 1 (y{E p ) cos(ntut)) + V^fy^Ep) cos 2 (nut)) E n 



V' 1 

p 



Vi (cosjlut) cos(nut)) + E n ( V{ {cos 2 (nut) cosjlut)) + ]-V Q ' (cos 2 (nut)) 



.1=1 



= V n /(2V P ) + %jV v + VijV p E n /A = A harm (E u ) + A inc (E aJ ,E n ) + A coh (E u ,E n ). 



(3) 



Using Esaki-Tsu characteristic, its derivative 

2V P 1-E 2 - E 
-. hi = — 

E r 



V'(E) 



(4) 



Ec (1 + E 2 Y 

and the definition of Fourier coefficients (0) , we calculate 



A h 



2V V 



(-1) 



1 



2k+l 



bE; 



2k+l 



\l + E z 



,1/2 



(5) 

(2k + 1 = n, k = 0,1,2,..., n = 3,5,7...; E u , n = 
E Ut n/E c ) for the harmonic component of total absorp- 
tion (see Fig. 0] upper), 



. 1 , 



En fW 



4K 



V' 



E„ 



,3/2 



for the incoherent component of total absorption, and 



(6) 



A. 



coh 



E„ 

4K 



In 



E n f 27r (l - El cos 2 x) cos(2nx) 
2 *" Jo (l + E 2 cos 2 x) 2 



dx 



(7) 



for the coherent component of total absorption. 

Note that Ai nc /E n , Eq. |@, is independent on n; it 
is always positive. The integral Q can be calculated 
analytically in two limiting cases E u <§C 1 and E u 3> 1. 
For instance, for n = 3 and E^ < 1, we expand V'(E) 
up to the term of E% and get A coh w -(7/64)^£ , 3 . In 
the limit £j w > 1 it easy to obtain A coh » -(2nE n )/ 'E% 
for any n. In general case the integral Q can be easily 
calculated numerically, see Fig. 0] lower. Note that A co h 
is always negative. 

It is also instructive to compare Ai nc and Aqq^ (se6 
Fig. . For i? w < 1 we have |A co/l /A mc | ~ £ 2 ™ < 1. 
In the opposite limit i? w S> 1, we have \A co h/ 'Ai nc \ ~ 
3> 1. Thus, a contribution of the coherent component 
to the total absorption is larger than a contribution of 
the incoherent component for a large pump amplitude. 
Of course, a contribution of both the coherent and the 
incoherent components is still less than a contribution of 
the harmonic component. 
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FIG. 4: (Color) Small-signal harmonic Ah [Eq. ©] (up- 
per subplot) and coherent A co h/E n [Eq. (J7J] (lower subplot) 
components of the absorption as a function of the pump field 
amplitude E^jEc and for different harmonics n = 3, 5, 7. 
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FIG. 5: (Color) Comparison of the incoherent and coher- 
ent absorption components. Dependencies of Ai nc /E n (black 
curve) and —A co h/E n (for different n) on the pump field am- 
plitude E w /E c . 
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APPENDIX B. DYNAMICS OF THE FIELD IN 
THE CAVITY 

Here we present the equations describing dynamics of 
electric field in the cavity and derive formula for the SL 
oscillator efficiency with an account of finite Q-factor of 
the cavity. 

For simplicity we consider a single-mode cavity with 
the eigenfrequency lo c tuned to the nth harmonic: u c ~ 
Lj n - We search the field in the cavity Efi e id in the form 

E fl eid = E n (t) cos[w„t + </> n (t)], (8) 

where the amplitude E n {t) and the phase <j>(t) are slow- 
varying in comparison to the carrier frequency u> n , i.e. 



with uj 2 pl = 4ire 2 N/(em ), m = (2fi 2 )/(Aa 2 ) and 
E n (t) = E n {t)/E c . Stationary solution of Eq. $H\i. E^, 
has the form 

E st UJ 2 ,T 

%- = —*-A{E ut E«), (14) 

where Q = u) n / r y c . This formula has been used to plot 
Fig. 

As follows from , the stationary value of the field 
inside the cavity is reached during the characteristic time 
~ 7c i.e. \E n /E n \ ~ 7 C . Therefore the condition of 
applicability of SVEA JjJJ is satisfied in a cavity with 
enough high Q: j c /u n = 1/Q < 1. 



\E n \ < U> n \E n \, \tj> n \ < Un\4>r, 



(9) 



Substituting (jHJ in the Maxwell equations and using the 
so-called slow varying envelope approximation (SVEA), 
which is valid in conditions ©, one can obtain 



dE n 
dt 



-E n 



2vr 



2tt 



(10) 



(11) 



Jc. s (k n ,t) = — [ dxR(x,k n )j c . s , (12) 
L Jo 

where j c is the phenomenological loss constant, R(x) is 
the cavity mode with wavenumber k n (R(x) ~ exp(ik n x) 
and u! c = ck n for travelling wave mode; R{x) ~ sm{k n x) 
for standing wave mode). In what follows we will ne- 
glect spatial dependencies of J c s because wavelength of 
microwave or THz radiation is larger than nanostructure 
length, k n L ~ L/X <C 1. Involved in Eq. I|12l) cosine and 
sine Fourier transformations of the current are defined as 

j c = 2j (v(t) cos uj n t), j s = 2j (v(t) sinw n i), 

where jo — eV p N and v(t) = V(t)/V p is the scaled ve- 
locity of an electron in SL arising under the action of ac 
field E{t) =E p {t)+E md {t). 

Next, we suppose that </> n (0) = and the exact reso- 
nance, uj n = lu c . Interaction of miniband electrons with 
the field E{t) is quasistatic and it is determined by the 
Esaki-Tsu dependence V(t) = y[i?(t)]. In these condi- 
tions we have j s (t) = j s (0) = 0, <j> n (t) = 4> n {0) = and 
j c = 2jqA. This approach means that for the quasistatic 
interactions of electrons with fields we can neglect disper- 
sion effects in the cavity. Now we have single equation 
for the field in the cavity 



dE n 
dt 



7cT 



E n 
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